Factorization of quantum charge transport for non-interacting fermions 
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We show that the statistics of the charge transfer of non-interacting fermions through a two-lead 
contact is generalized binomial, at any temperature and for any form of the scattering matrix: 
an arbitrary charge-transfer process can be decomposed into independent single-particle events. 
This result generalizes previous studies of adiabatic pumping at zero temperature and of transport 
induced by bias voltage. 
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I. INTRODUCTION 



Noise is usually considered as an unwanted characteristics of electronic circuits. Sometimes, however, studying 
electronic noise gives important information about physical systems. In particular, noise in small junctions is affected 
by the fermionic statistics of electrons and is an essential part of quantum mesoscopic transport (see, e.g., Ref. fT|. 
Even more information about the quantum behavior of charge carriers is encoded in the full probability distribution 
of the transferred charge: a concept known as full counting statistic^. 
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The full counting statistics of electrons is most studied for the simplest model of non-interacting fermions in various 
setups. In this model, the general expression for the proba bility distribution of the transferred charge is given by 
an elegant determinant formula of Levitov and LesovikpEEEl This formula expresses the generating function for the 
probability distribution as a certain functional determinant involving the single-particle scattering matrix. An exact 
calculation of this functional determinant may be performed in several special setups of the charge transfer^'^ '^'^*"'^, 
and until recently most of the results on the full counting statistics in this model concentrated on studying such 
exactly solvable cases. 

However in the recent years a progress has been made in understanding general properties of charge transfer encoded 
in the determinant formula. Namely, it has been shown that, in the case of a contact with two external leads, the 
total electronic transfer is given by a superposition of uncorrelated elementary charge transfers of single electrons. 
In each of such single-electron events, one electron has a certain probability p to pass through the contact (thereby 
transferring one quantum of charge) and the probability 1 — p to reflect (resulting in no charge transfer). This 
decomposition into elementary charge-transfer events (dubbed "generalized binomial" statistics^ has been shown 
under various assumptions: at zero temperature for a charge transfer driven by a time-dependent bias voltag^, at 
zero temperature for the adiabatic-pumping problerrPSl and in the wave-packet formalism (a finite number of wave 
packets with a momentum-dependent scattering^. 

The goal of the present paper is to lift the assumptions made in those previous works and to show that the 
generalized binomial statistics is universally valid for any charge-transfer problem involving non-interacting fermions: 
at any temperature and for any time- and energy-dependent scatterer. A necessary condition remains that the initial 
state does not involve any entanglement between the leads or between different charge sectors in each lead. 

It may be therefore practical to describe the statistics for noninteracting fermions by the distribution /i(p) of 
the effective transparencies of the elementary charge transfers, instead of the probabilities of the total transfer. To 
illustrate this construction, we list several known examples where the distribution fi(p) is exactly known. Finally, to 
illustrate the importance of the non-interacting assumption, we also present several examples of interacting systems 
where the full counting statistics is not generalized binomial. In those cases, we may interpret the effect of interaction 
as a shift of the transparencies p to the complex plane. 

The paper is organized as follows. First, to fix the notation, we briefiy review the derivation of the Levitov-Lesovik 
determinant formula in Section |ll] In Section III we obtain our main result: a decomposition of an arbitrary charge 
transfer process into independent tunneling events. In Sections |IV| and \V\ we illustrate our construction with several 
noninteracting and interacting examples, respectively. Finally, we conclude with Section [yT) where we discuss possible 
generalizations and applications of our results. Some technical details are relegated to Appendices A, B, C. 



II. GENERATING FUNCTION AND DETERMINANT FORMULA 



In this section we briefiy review some known results from the theory of full counting statistics for non-interacting 
fermions (mainly following the method of Ref. 11). We also introduce the notation and specify the assumptions used 
in the subsequent sections. 



A. Generating function for full counting statistics 

A typical full-counting-statistics problem involves determining the probabilities of a given charge transfer through 
the junction. In this paper, we will consider the case of a junction with two leads, and therefore the definitions in 
this section will be written for a two-lead junction, even though they admit a natural generalization to a many-lead 
situatiorf^. In the case of two leads, the charge transfer is characterized by one number q: an integer number of 
electrons moved from the left to the right lead (or —q electrons from the right to the left lead). The ingredients of the 
problem are: (i) the initial density matrix of the system po and (ii) the unitary evolution operator U (over the time 
when the measurement is performed). To formally define the transferred charge q, we also need (iii) an operator Q 
which counts the particles in one of the leads (say, the right lead). 

Then one easily sees that the probabilities of a given charge transfer q can be found from the generating function 

X(A) = Tr (^poU^e^'QUe-^^Q) / Tr po , (1) 

where A is an auxiliary variable and the trace is taken in the many-body Fock space (for technical reasons, it is 
convenient to explicitly include the normalization of the density matrix). For a good definition of charge transfer 
probabilities, one needs to further assume that different charge sectors are not entangled, i.e., that pQ commutes with 
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Q. Under this assumption one finds that thus defined x(A) may indeed be interpreted as a generating function 



X(A) = E 



g— — oo 

for some probabilities < < 1. The periodicity of x{^) (or, equivalently, the quantization of q) follows from the 
fact that the spectrum of Q is integer. It will therefore be convenient for our purposes to treat x(A) as a function of 
the complex variable 

u = (3) 

Then Q becomes the Laurent series of the complex function The Fourier components of x(A) (the coefficients of 

the Laurent series) define the probabilities Pq, while the derivatives at A = give the moments of the charge transfer 
('?")• 

In the case of a periodic pumping, it may be useful to define a "full counting statistics per period" ^ as the extensive 
part of xW'- 



Xo(A) = exp 



lim 



1 



Inx(A) 



p 



(4) 



where the full counting statistics x('^) is collected over Np periods. Defined in this way, the characteristic function 
Xo(A) may be used for calculating extensive quantities (such as cumulants of the transferred charge), but the effects 
associated with the beginning and ending of the measurement!^ are neglected. 



B. Determinant formula for noninteracting fermions 

In the case of non- interacting fermions, the general expression ([T|) may be sim plified and rewritten in terms of 
single-particle operators: the so-called Levitov-Lesovik determinant formulaP^'^'^. In this section we briefiy present 
its derivation following the approach of Ref . lit 

For non-interacting fermions, the multi-particle operators (acting in the Fock space) are related to the corresponding 
single-particle operators (acting in the Hilbert space of one particle) with the help of the second-quantized formalism. 
For notational convenience, we will use the "hatted" notation (as in Eq. [T|) for operators acting in the Fock space, while 
the "unhatted" letters will denote the corresponding single-particle operators. For quantum-mechanical observables 
(such as the charge operator), the correspondence is given hy Q = tp^Qip (where '0^ and V' are the fermionic creation 
and annihilation operators, and the summation over single-particle states is assumed), while for their exponentials 
(such as the evolution operator or the density matrix) the correspondence is IJ — exp[?/;^ (In [/)■)/)] (see Appendix [a| for 
details) . 

Now we are ready to derive the Levitov-Lesovik determinant formula, starting from ([T]) and making the following 
assumptions: 

1 ) The operator Q is a fermionic bilinear. This is obviously true for the operator of the charge in the right lead. 
However, more generally, one can use the same formalism for studying full counting statistics for other physical 
quantities (e.g., spin). 

2) The evolution operator U is an exponential of a fermionic bilinear. This is equivalent to requiring that the 
Hamiltonian is a fermionic bilinear (see Appendix [A| , i.e., that the fermions are non-interacting. 

3) The density matrix po is an exponential of a fermionic bilinear, po — e~''°, where Hq — ip^hotp is a Hermitian 
fermionic bilinear (note that, if using Eq. ([T]), we do not need to normalize the density matrix). Usually, one considers a 
thermal equilibrium with Hq given by the physical Hamiltonian divided by the temperature. However, our assumption 
is much less restrictive. For example, one can take an initial state with different temperatures of the leads and/or of 
different channels or even prepare an equilibrium state using a different Hamiltonian. 

4 ) The operators Q and po commute. This condition expresses the absence of charge entanglement in the initial state 
and is a necessary condition for interpreting x(A) as a generating function for charge-transfer probabilities. 

The trace of the product of exponentials of fermionic bilinears in ([ij can be rewritten as a determinant of an 
operator acting in the single-particle Hilbert space (see Ref. ITT] and Appendix [A| for details): 

X(A) = det [l + npiU^e'^'^Ue-'^'^ - 1)1 . (5) 
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Here L'^ is a single-particle evolution operator and 



is the occupation- number operator. 



IT'F = 7 = J— (d) 



C. Regular ization of the determinant 



Depending on the dimensionality and on the asymptotic behavior of the operators involved in (|5| , the determinant 
may require a regularization. 

If the dimension of the single-particle Hilbert space is finite (this is the case in the wave-packet formalism of Ref . |S| , 
all the matrices in ([s]) are finite-dimensional, and no regularization is needed. 

If the single-particle Hilbert space has an infinite dimension, then a regularization may be needed or not, depending 
on the properties of the Hamiltonians describing the evolution (in the operator U) and the initial state (in the operator 
p). Usually, one considers a certain dispersion of the particle in the leads (the same in the evolution and the initial- 
state Hamiltonians), with the fermions filling the states up to some Fermi energy. The determinant ^ is well-defined, 
if the matrix tends to unity sufficiently fast at both large positive and negative energies. At large positive energies, 
np tends to zero, and the matrix under determinant tends to unity. At large negative energies, the determinant is cut 
off either by the Fermi energy or by the energy dependence of U: at very low energy, the particle cannot penetrate 
through the barrier, then U commutes with Q, and the matrix again tends to unity. So the determinant ^ is well 
defined if one uses exact dispersion relation and the energy dependence of the scattering matrix. 

The only (but historically very common in literature) case when the determinant needs a regularization is the ap- 
proximation where the spectrum is linearized (so that the Fermi energy does not serve as a regularizing paramet er) and 
the scattering matrix is assumed to be energy-independent ("instant scattering" or "adiabatic pumping" Ijniit)|21i21i4| 
In that case, the matrix in (|5| is infinite-dimensional and tends to a non-unity matrix in the negative-energy asymp- 
totics, therefore the determinant needs to be regularized. The regularization may be performed by admitting a weak 
energy dependence of U, so that the determinant converges, and then letting the corresponding energy scale tend 
to infinity. Technically, this procedure amounts to simply re-expressing the determinant in a manifestly convergent 
form where the weak energy dependence of U becomes inessential. This approach was used in Ref. If and then more 
generally (but also more formally) in Ref. il5t For completeness, we review the details of adiabatic approximation in 
Appendix |B] 

Therefore, we assume, without loss of generality, that an appropriate energy dependence is already included in the 
evolution operator U, and the determinant ([5| does not need any further regularization. Moreover, in the future 
discussion we will be interested in the zeros of the determinant, which are regularization independent. 



III. DECOMPOSITION INTO SINGLE-PARTICLE EVENTS 



A. Zeros of the generating function 

To demonstrate that the full counting statistics for non-interacting fermions can be decomposed into elementary 
single-particle transfer events, we first re-express the determinant (|5| in terms of the spectral properties of a single 
operator, whose eigenvalues give transmission probabilities of those elementary transfers. 

This transformation proceeds as follows. Using the fact that the charge operator Q is a projector (Q^ = Q), we 
rewrite e*'^'^ = 14- (e*'^ — 1)Q in ([s]), in order to obtain, after some simple algebra, 

x{X)^det{[l + {e^^-l)X]e-^^Q) . (7) 

Here we have introduced a new operator, 

X={l-nF)Q + nFU'fQU. (8) 
One can see that the spectrum of the operator X is 
1. Real, 



2. Confined to the interval [0, 1]. 
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Both properties become obvious if one considers a Hermitian operator conjugate to X by a similarity transformation: 

X = {nFy^^^Xinpy/^ = (1 - nF)Q + {npf'^U^QUinFY'^ (9) 

(we have used our assumption [np^Q] = 0; np is Hermitian and positive definite so that [np)^^^'^ is well definecP^. 

The reality of the spectrum follows from the hermiticity of X. The constraint on the eigenvalues follows from the 
observation that both X and \ — X are represented as sums of two positive operators (see Eq. ^ and the one obtained 
from it by replacing Q I — Q). 

The spectral properties of the operator X immediately imply a constraint on the positions of zeros of the generating 
function x(A) analytically continued to the complex plane of the variable u = e*^, 

X{u) =det{[l + {u-l)X]u-^) . (10) 

The zeros of x(m) can be 0, oo (from the factor u~'^) or coincide with 1 — X^^, where X„ are the eigenvalues of the 
operator X. From the spectral properties of X, it then follows that zeros of x(w) are confined to the negative real 
axis of the complex u-plane. This conclusion has been previously reached in Ref. ITU] under the assumptions of zero 
temperature and instant scattering. Generalizing it to arbitrary temperatures and to arbitrary evolution operators 
(under the assumptions formulated in Section constitutes the main result of our present paper. 



B. Generalized binomial statistics and effective transparencies 



The expression ( |10[ ) may be further interpreted in terms of decomposing the charge transfer into a superposition 
of independent tunneling events, with the eigenvalues of X giving the probabilities of those tunneling processes. To 
make this interpretation more explicit, we get rid of the last factor exp[—iXQ] in ([To| bjEH 



X(A) = det([l + (e-^^-l)X]e^*^^e*^^e-*^'^) 
= e'^t'^^^-^Met ([1 + (e^^ - 1)X] g-*^^) 



(11) 



[1 + (e'^ - 1)X„] , 



where in the last line the product is taken over the eigenvalues of the operator X, and 

{q) = -tdxlnxiX) = ti- (X - Q) = ti- [np{U^QU-Q)] 
\=o 



(12) 



is the average total transferred charg^i^ 

Each factor [1 + (e*'^ — l)Xn] in (111 corresponds to an elementary single-electron event with the transmission 
probability X„ and reflection probability 1 — Xn- The other exponential factors in (111 correspond to a deterministic 
"background" charge transfer, which does not produce any noise. 

Such a statistics given by a superposition of individual transmission events with different transparencies was dubbed 
generalized binomial statistics in Ref. 8 Note that the effective transparencies are given by the spectrum of the 
effective-transparency operator X and depend in a complicated way on the time dynamics of the junction and on the 
initial thermodynamic state.'^ 



C. Spectrum of the effective-transparency operator 

Using the result derived above, we may now describe the full counting statistics of non-interacting fermions by the 
distribution of effective transparencies (the density of eigenvalues of the operator X), instead of the probabilities Pq. 
Following Ref. (THl we define the effective-transparency density as 

M(p) = tr<5(p~X). (13) 

Since the spectrum of X defines the density of zeros of x(") oii the negative real semi-axis of u-plane, the density 
of states ii{p) may be related to the jump of the derivative of xW at the corresponding point. A simple calculation 
gived^ 

Jt-l-iO I 

= — Im dp In x{p — iO) , (14) 



2ttv(1 — v) 



dxlnxW 
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where 

P-T^-T^- (15) 

Conversely, knowing the spectral density one can reconstruct the cumulant generating function Inx(A) (up to 
an overall deterministic charge transfer) by using Q: 

lnx(A) = -iATrQ+ / In [l + (e*^ - l)p] . (16) 

Jo 

As for the structure of the spectrum of effective transparencies, there are two alternatives: the spectrum may be 
either discrete or continuous. An analysis of available examples gives the following possibilities: 

1. In a finite-dimensional problem (e.g., a finite number of incident wave packet Jsl)^ the matrix is finite-dimensional, 
and the spectrum of effective transparencies is obviously finite and discrete. 

2. At zero temperature, in the instant-scattering approximation, for a continuous and periodic time dependence of 
the scattering matrix^*^, the spectrum is discrete, with possible accumulation points ai p — and p — 1. 

3. At zero temperature, in the instant-scattering approximation, for a discontinuous time dependence of the scat- 
tering matrix, the spectrum typically becomes continuou^iSl 

4. From the above two examples, one can deduce that at zero temperature, in the instant scattering approximation, 
for a continuous non-periodic time dependence (a finite-time pumping pulse), the spectrum is discrete, if the 
scattering matrix returns to its initial value, and continuous otherwise. This can be shown by a linear fractional 
mapping of the time axis on a circl^, thus reducing the non-periodic problem to a periodic one, with either 
continuous or discontinuous time dependence of the scattering matrix. 

5. For finite temperatures or energy-dependent scattering, the spectrum of X is typically continuous. This can be 
seen in the examples presented in Section |IV[ The physical reason for a continuous spectrum is that particles 
incident at different energies have different effective transparencies, and thus integration over energies smears 
the spectrum. 

If the spectrum fi(p) is continuous, the singularities form a branch cut of the function x('^) along the negative real 



axis (or along a part of it). At this branch cut, the function x(m) is discontinuous, with the jump obeying (14 1. A 
more detailed analysis of the spectrum of the effective transparencies will be presented in a subsequent public at ioiP^. 

Note that in the case of zero temperature and adiabatic pumping, the spectrum of effective transparencies may 
alternatively be found from a different operator labeled rip in Ref. 10 and M in Ref. 1^ It is worth emphasizing 
that although the two operators are different, at zero temperature and in the instant-scattering approximation their 
spectra coincide, as shown in Appendix [C| 



IV. EXAMPLES: NON-INTERACTING FERMIONS 



In this section, we illustrate our discussion with several known examples. 



A. T — 0. Optimal (Lorentzian) pulses 

An instructive example where the full counting statistics is exactly computable is the case of several Lorentzian 
voltage pulses of one-quantum intensity each: J V{t) dt = ±27rS. In this case, as shown in Ref. [SJ the full counting 
statistics involves only a finite number of transferred electrons. For the effective-tr ansp arency spectrum, this implies 
a finite discrete spectrum with the number of levels equal to the number of pulseaSES. In the special case when all 
pulses are of the same polarity (this also includes the setup of a constant voltag^ as a limiting case), all the effective 
transparencies are degenerate and equal to the transparency of the junction g: 

^^{p)^WS{p-g), (17) 

where W = (27rfi.)~^ J V{t) dt is the number of pulses. 



7 



0^ 
0.0 

















p 



0.2 



0.4 



0.6 



0.8 



1.0 



FIG. 1: The effective-transparency distribution (18 1 for g = 1 (blue dashed line) and g = 0.96 (red solid line), plotted in the 
units of 4G/7r^. For g = 0.96, there is a spectral gap between p = 0.4 and p = 0.6. 



B. T = 0. Periodic opening of a quantum point contact 

An interesting exactly solvable example with a continuous spectrum of effective transparencies has been consid- 
ered in Ref. |TH] a periodic train of openings and closings of the contact. In this case, the distribution of effective 
transparencies is (Fig. [T]) 



G 



Re 



2p\ 



7r'p{l-p) V(l-2p)2-(l-5)' 



where the overall intensity [per opening/closing cycle, as in (H|] is 



G = 21n 



TTWo 



(18) 



(19) 



Here w is the fraction of a period during which the contact is open (with the transparency g) and wq is a regularization 
given, essentially, by the fraction of a period taken by the switching time. 

One can see that if the contact is not fully open {g < 1), then there is a spectral gap around /i(p) = 1/2. We will 
comment on this feature in a future public at iorP^l. 

C. T 7^ 0. Zero bias voltage 

Remarkably, the same formula ([T8| gives the distribution of effective transparencies in the case of a finite temper- 
ature T (and zero bias voltage)** The overall intensity is then given by 



G = 



■KtT 



(20) 



and the observation time t is assumed to be much longer than fi/T, so that only the extensive part (0 is retained. 



EXAMPLES: INTERACTING FERMIONS 



The derivation of Section [III] shows that the absence of interactions between fermions leads to a generalized binomial 
statistics of the charge transfer between two conductors. In this section, we show that if the fermions interact, the 
zeros of the generating function x(w) may shift to the complex plane. 

Before presenting examples of full counting statistics for interacting fermions, let us use a two-particle example to 
gain some intuition about the positions of the roots of xi^)- Consider a system of two particles which propagate 
from the left to the right lead of the contact. The full counting statistics is then given by a quadratic polynomial 
x{u) = ^0 + Piu + P2U^ with the probabilities Po,i,2 obeying Pq + ^'i + ^2 = 1- The roots of the polynomial become 
non-real if 



Pi < 4P0-P2 , or, equivalently, 



^P2 > 1 



(21) 
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i.e., if the probability Pi of separating the two particles becomes too small. If one considers a strong attraction 
between particles, so that they become inseparable, then Pi = 0, and the roots of xC"*^) become purely imaginary. One 
thus sees that in this example (i) the non-reality of roots may be loosely related to an attraction between particles; 
(ii) a certain threshold of interaction may be needed to shift the roots of x('") from the real axis. 

We now turn to several more complicated examples of interacting systems where the full counting statistics may 
be exactly computed. 



A. Charge transfer in a normal metal — superconductor point contact 

A more general example of an electronic system with attraction leading to charge transfer quantized in pairs of 
electrons (and hence to a shift of the zeros of xi''^) from the real axis) is considered in Ref. [151 a contact between a 
normal metal and a superconductor. In the low-temperature limit and at low bias voltage V (so that T ^ eV <C A), 
single-particle tunneling is suppressed, and the full counting statistics is given by 

x{u)^{l-gA + 9Au')'^ . (22) 

Here gA is the probability of Andreev reflection and W = eVt/{2Trh) is the effective number of transfer attempts 
during the observation time t. Similarly to our two-particle example above, the zeros of x('") Eire in this case purely 
imaginary. 



B. Two particles scattering on a resonant quantum dot 



Another example of an interacting two-particle scattering problem has been studied in Ref. 20 In that work, one 
considers a quantum dot characterized by its resonances and by a Coulomb interaction between the particles. In this 
model, the scattering matrix is exactly computed and then used to find the full counting statistics in a two-particle 
example. For the example considered in that work (a singlet two-particle wave packet in the shape of an exponentially 
truncated plane wave), the probabilities to transfer one and two electrons are given by (we refer the reader to Ref. [20l 
for details of the model and the derivation) 



Pi = 

P2 = 



2(3 



(l + 3/3)(l + /3)2 
1 

(l + 3/3)(l + /3)2 



2 + 3/3 
1 



1 



(52 + l_ 

s/3 [452 + 3(3 -1-/3)] 



(23) 



2{l + I3){S^ + l){d^ + s^) 



where /? is the dimensionless dwelling time on the dot characterizing the resonance, S characterizes the strength of the 
Coulomb interaction (it is related to the parameter a introduced in Ref.EOlas d = a/3/{l+f3)), and s = (3-|-/3)/[2(l-|-/3)]. 

Using the criterion ( |2ip , we can find the range of parameters S and /? where the statistics is generalized binomial 
(see Fig. [2]). Depending on the dwelling time, an interaction may either shift the roots of x('^) from the real axis or 
not. 

Although this is a two-particle example, a qualitatively similar behavior of roots can also be found in multiparticle 
systems. For example, in Refs. 1211221 full counting statistics has been calculated for charge and spin transport through 
a Kondo dot. In that problem, one can have both regimes with roots real and complex, depending on the parameters. 



C. Charge pumping through a single-electron transistor 

An interesting example with interaction where the statistics is nevertheless generalized binomial is studied in Ref. H31 
The authors of that work consider a charge pump based on a nearly open quantum dot in the Coulomb blockade 
regime. In spite of the interaction on the dot, the problem is mapped onto an effective quadratic fermionic Hamiltonian 
(including a Majorana fermion representing the resonant level at the dot). Therefore, our formal conclusion about 
the generalized binomial statistics remains valid for this interacting system as well. 

Indeed, one can easily verify that the generating function (per period) found in Ref. 1231 

lnx(A) = -.A-H^/&ln(l--, 



n_(l — n+) n_|-(l— n_) 
P 1 



(24) 
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generalized binomial statistics 






complex roots 


6 



12 3 4 



FIG. 2: The shaded area shows the region in the parameter space of 5 and (5, where the roots of x('m) are complex and the 
statistics is not generalized binomial. 
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FIG. 3: The effective-transparency distribution (251, in the units of F/hui (thick red line). vanishes abruptly at p 

{hujf/({hujf + r^) (shown for huj = F). 



contains only real roots p (at any temperature). Here w is the pumping frequency, n±{e) — (e^'^^'^'''^^/^ + are 
the Fermi occupation numbers at the temperature T, T is an energy scale characterizing the strength of the Coulomb 
interaction in the dot, and p = (1 — e*^)~^, as usual. 

In the zero-temperature limit, the effective-transparency distribution ^{p) can be easily calculated analytically: 

^(^)= r for ^''^ 

I" lor p > (huiY+T^ ■ 

The spectrum is continuous with a gap at s^ ji^^y^ < p < 1, as shown in Fig. jsj (one can check that the gap survives 
at finite temperatures, but changes its magnitude and shifts towards the center of the interval). 



VI. CONCLUSION 



To summarize, we have shown that the charge-transfer statistics for a system of non-interacting fermions is gener- 
alized binomial: it factorizes into individual single-particle tunneling events. This property holds under very general 
conditions: at arbitrary temperature and for an arbitrary time- and energy-dependent scatterer. Note that the con- 
verse is not true: an interacting system may either obey or disobey this factorization property, as we have seen in 
examples. 

The factorization property may be formulated in terms of the positions of zeros of the generating function x(u) 
(or, equivalently, of the singularities of Inx(u)). The factorization implies that those singularities are confined to 
the negative real axis. In this case one can fully characterize full counting statistics by the distribution of effective 
transparencies /i(p}. All charge-transfer properties (moments and cumulants of the transferred charge, entanglement 
between the leads^^, etc.) can be expressed in terms of /i(p). 

The absence of factorization corresponds to a shift of singularities from the real axis. In this case, x(m) can still 
be characterized by its zeros albeit in the complex w- plane. This description is somewhat similar to the approach of 
Lee and Yang in statistical physics, where a partition function is characterized by the zeros of its analytic extension 
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to complex values of physical parameters such as magnetic field or chemical potential^. In statistical physics, this 
approach turned out to be very useful for understanding possible phase transitions: a phase transition is determined 
by the location of zeros of the partition function relative to the real (physical) axis of parameterd^. From this 
perspective, there remain interesting questions of characterizing interactions in the scattering problem by the positions 
of singularities of Inx(u) and of understanding possible physical implications of their shift away from the real axis. 

As an interesting extension of our results, one can consider constraints on full counting statistics for non-interacting 
electrons in a multi-terminal setup. In this case, a certain "reality condition" on the zeros of the generating function 
is possible^*. However, the factorization property does not hold, being specific to the two-terminal case. 



Another possible generalization of our construction is to use the formalism developed in Sections |II| and III for other 
types of problems with quadratic Hamiltonians, not necessarily involving physical charge transfer. As an example of 
such an application, one can study the statistics of staggered magnetization of a XY spin chain (mappable to free 
fermions by a Jordan- Wigner transformation). See also Ref. 'SB'for applications of those ideas to analyzing the entropy 
in such systems. 
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APPENDIX A: FERMIONIC BILINEARS AND TRACE IDENTITIES 

For reader's convenience, we present in this appendix some helpful identities on the exponentials of fermionic 
bilinears. Those identities are then used to derive the determinant formula A more detailed account including 
proofs may be found in Ref. ITTI 

For any single-particle observable A, the corresponding multi-particle operator acting in the Fock space is given by 
the fermionic bilinear 

A = ijj'^AiP . (Al) 

This construction provides a representation of the gl{n) algebra in the 2"-dimensional Fock space (n denotes here the 
dimension of the single-particle Hilbert space). Indeed, one easily checks that commutators of fermionic bilinears are 
given by 

[A,B]^[A^]. (A2) 
Exponentiating these relations, one obtains a representation of the group GL{n)^^ 

e^e^^e^, if e^e^ = e^. (A3) 

One can use this representation to show that the evolution operator for a system of noninteracting fermions is 
an exponential of a fermionic bilinear. Indeed, the multi-particle evolution operator is given by the time-ordered 
exponential 



U = Texp 



-i j H{t)dt 



(A4) 



where the time-dependent Hamiltonian is a fermionic bilinear. This time-ordered exponential can be represented 



as a chronologically ordered product of (infinitesimal) exponents of fermionic bilinears and, by virtue of (A3 1 may be 
rewritten as a single exponential 

[/ = exp[?/'^(lnC/)V'], (A5) 

where 



U = Texp 



-i j H{t)dt (A6) 
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is the single-particle evolution operator. 

For traces of such exponentials one has a simple formulaP^ 

Tre^ = det(l + e-^), 



(A7) 



where the trace is taken in the Fock space, and the determinant is in the single-particle space. For traces of products 
of exponentials of fermionic bilinears, one finds, using (A3) and (A7l: 



Tr 



det(l 



Ai A2 



From this relation, a derivation of the Levitov-Lesovik determinant formula easily follows: 



Tr 



-iXQ 



Trpo 

where n^? is defined by (|6l 



dot (1 + e-''"We'^^Ue-'^Q) 
det {1 + e-f^°) 



= det (1 + np {U^e'^^Ue-'^Q - l)) 



(A8) 



(A9) 



APPENDIX B: ADIABATIC AND INSTANT-SCATTERING APPROXIMATIONS 



The determinant formula ([5| has been derived under the very general assumptions specified in Section [IT] A phys- 
ically relevant and technically simpler approximation is possible in the situation when the two leads involve asymp- 
totically free electrons, so that the scattering may be described in terms of an instantaneous scattering matrix S{t, e). 
This matrix relates the asymptotic states for electrons at energy e scattering on the instantaneous Hamiltonian H{t). 
If the Hamiltonian varies slowly on the scale of the "scattering time" hdlogS/de, then the evolution operator U 
may be approximated by the Wigner transformation of S{t,e)^. The corresponding adiabaticity condition may be 
formulated as 



dt de 



< 1. 



(Bl) 



At zero temperature, this condition allows to neglect the energy dependence of S{t,e), and to replace the scattering 
operator in ([5]) by the instantaneous scattering matrix S{t) taken at the Fermi energy e (since only the states in 
the vicinity of the Fermi energy contribute to the determinant ^P^'^. 

At finite temperature, to neglect the energy dependence of the scattering matrix, one also needs that the energy 
dependence of S is small on the scale determined by the temperature. In addition to the adiabaticity condition (Bl I 
we require: 



T 



S 



_i dS 



< 1. 



(B2) 



The two conditions (Bll and ( |B2| constitute the instant- scattering approximation (equivalent to the adiabatic ap- 
proximation in the case of zero temperature). In this approximation, the determinant formula becomes 



X(A) = det (1 + npiS^ e'^^ Se-'^Q - 1)) 



(B3) 



where S{t) is a local in time unitary scattering matrix. 

The determinant (B3l, in turn, requires a regularization at negative energies, since the operator under the determi- 
nant tends at infinite negative energies to S^e^^'^ Se^^^^ , a non-unity matrix. This regularization may be performed 
either by taking into account the finiteness of the negative-energy spectrum or by re-introducing a weak energy de- 
pendence of S{tj e) so that lim^^-oo S{t, e) = 1. In the latter case, the easiest way to implement the regularization is 



to rewrite ^ 



aa^ 



X(A) = det ( [1 + npiU^'e'^'^Ue-'^Q - 1)] • e'^"''0c/te-'^"^'3c/) 



(B4) 



This manipulation is obviously admissible for the "full" evolution operator U corresponding to the regularized scat- 
tering matrix with lmig^_oo S(t,e) — 1. However, in the new expression (B4|, one can simply replace [/ by a local in 
time scattering matrix S(t) to obtain a fully convergent expression in the instant-scattering limit. 

The regularized formula ( |B4[ ) proposed in Ref.[15]is equivalent to the previously existing regularization prescriptions. 
Our discussion in this paper does not rely on this regularization, since we always consider the general case of an 
arbitrary non-interacting evolution operator and therefore may assume that the negative-energy asymptotics is suitably 
regularized. 
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APPENDIX C: SPECTRUM OF THE EFFECTIVE TRANSPARENCIES IN THE ADIABATIC LIMIT AT 

ZERO TEMPERATURE 

Generalized binomial statistics for non-interacting fermions has already been proven in our earlier publication under 
the assumption of zero temperature and instant scattering Of course, that result now follows from the more general 
argument of the present work. However, we find it instructive to demonstrate explicitly the equivalence of the two 
results under those more restrictive assumptions. 

It has been shown in Ref. 10 that, at zero temperature and in the instant-scattering limit, the effective transparencies 
of elementary single-particle events are given by the spectrum of the operator 

n'p = QSnpS^Q (CI) 

(the same operator was denoted M in Ref. I16p . In this Appendix, we explicitly show that the spectra of n|n and of 
X coincide, except for the eigenvalues zero and one. 

Indeed, at zero temperature np is a. projector operator, and therefore the "hermitized" operator X given by 
whose spectrum coincides with the spectrum oi X, is block-diagonal. The block at positive energies equals Q and 
only gives eigenvalues zero and one. The block at negative energies can be written (at zero temperature) as 

X^ = npS^QSnp = R , where R = QSnp (C2) 



(we have also used that Q is a projector). On the other hand, one can rewrite (CI I as 



= RR^ , (C3) 

from where it becomes obvious that the spectra of X_ and rip coincide, except for zero eigenvalues (their eigenvectors 
can be related to each other by the operators R and i?'''). This completes the proof of the coincidence of the spectra of 
X and Up (modulo eigenvalues zero and one, which do not contribute to the noise, but only to the overall deterministic 
charge transfer). 
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